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基于 改进 自 适 应 黑洞 机 制 的 引力 搜索 算法 
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摘 要 :针对 基本 引力 搜索 算法 (gravity search algorithm,，GSA) 易 早 训 、 易 陷入 局 部 最 优 、 缺 少 有 效 加 速 机 制 等 缺点 ， 
提出 基于 改进 自 适 应 黑洞 机 制 的 GSA(Improved adaptive black hole gravity search algorithm，IABHGSA)。 通 过 改进 
Tent 映射 对 种 群 初始 化 ， 使 得 初始 种 群 的 分 布 更 随机 、 均 匀 、 人 遍历 ， 增 强 算 法 的 全 局 勘探 能 力 ; 引入 改进 自 适 应 黑 
洞 机 制 ， 根 据 粒子 进化 情况 选择 位 置 更 新 策略 ， 使 得 位 置 更 新 更 为 合理 ， 有 效 减 小 粒子 陷入 局 部 最 优 的 可 能 性 ; 通 
过 基于 学 习 思 想 的 最 优 与 最 差 粒 子 更 新 策略 ， 增 强 算法 逃离 局 部 最 优 的 能 力 以 及 提高 算法 的 寻 优 速度 ; 引入 群体 迁 
徙 ， 为 算法 提供 有 效 的 加 速 收 全 机制。 最 后 ， 选 取 8 个 基准 测试 鸥 数 对 IABHGSA 进行 测试 ， 并 与 相关 算法 的 实验 
结果 进行 对 比 ， 结 果 证 明 IABHGSA 有 更 好 的 寻 优 性 能 。 
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Gravity search algorithm based on improved adaptive black hole mechanism 


Xu Wenjun, Wang Xihuai', Xiao Jianmei, Gu Junyu 
(College of Logistics Engineering, Shanghai Maritime University, Shanghai 201306, China) 


Abstract: Aiming at the shortcomings of the basic Gravity search algorithm (GSA) , such as prone to premature maturity, 

easy to fall into local optimum, and lack an effective acceleration mechanism, this paper proposed a GSA (IABHGSA) based 

on the improved adaptive black hole mechanism. The algorithm used the improved Tent mapping to initialize the population, 

which makes the distribution of the initial population more random, uniform and traversal, and it enhanced the global 

exploration capability of the algorithm; The algorithm introduced an improved adaptive black hole mechanism and selected 

the position update strategy, according to the evolution of the particles, which makes the position update more reasonable, 

Effectively reduce the possibility of particles falling into local optima; Through the optimal and worst particle update strategy 

based on learning ideas, enhanced the algorithm's ability to escape from local optima and improved the algorithm's 

optimization speed; The algorithm introduced group migration to provide effective accelerated convergence mechanism for 

the algorithm. Finally, this paper selected 8 benchmark functions to test IABHGSA and compared with the experimental 

results of related algorithms, the results show that IABHGSA has better optimization performance. 

Key words: gravity search algorithm; 1mproved tent Mapping; adaptive strategy; particle position update; group migration 
0 ”引言 等 人 在 GSA 的 位 置 更 新 公式 中 引入 Levy 飞行 策略 ， 改进 后 
的 GSA 具有 更 好 的 寻 优 精度 和 稳定 性 ， 但 是 并 未 明显 提高 

引力 搜索 算法 (GSA) 是 Rashedittl 等 人 受到 物理 学 中 万 有 ”GSA 的 收敛 速度 ; MirjalilitSj 等 人 将 “向 全 局 最 优 移动 ”项 加 

引力 定律 的 启发 而 提出 的 智能 优化 算法 。 GSA 通过 粒子 间 引 ”入 GSA 的 速度 更 新 公式 ， 提 高 GSA 的 收 义 速 度 和 逃离 局 部 
力 的 相互 作用 来 交流 信息 ， 其 概念 简单 ， 容 易 实现 且 算法 收 ” 最 优 的 能 力 ， 然 而 GSA 可 能 因为 最 优 粒子 的 错误 引导 而 早 
敛 性 好 证 。 研 究 表明 ，GSA 的 寻 优 精度 和 收 和 敛 速 度 对 比 粒 子 ” 熟 收 勾 ; Su 等 人 在 GSA 的 引力 系数 中 加 入 基于 粒子 适应 度 
群 优 化 算法 (Particle Swarm Optimization，PSO) 、 遗 传 算 法 值 的 反馈 机 制 ,并 加 入 精英 记忆 保留 策略 ,提高 了 GSA 的 收 
(Genetic Algorithm, GA) 等 其 他 优化 算法 具有 一 定 的 优势 止 。 敛 速度 和 重 棒 性 ， 但 GSA 的 寻 优 精度 对 比 改进 前 提升 有 
因此 在 优化 调度 是 、 参 数 辨 识 外 、 函 数 优化 外 等 领域 , GSA 得 ” 张 维 平 09 等 人 在 GSA 中 引入 反 向 学 习 策略 、 精 英 策略 和 边 
到 广泛 的 应 用 。 界 变异 策略 ， 显 著 地 提高 了 GSA 的 全 局 探索 能 力 和 局 音 
GSA 与 常见 的 生物 启发 式 算法 相似 , 都 属于 群 智能 优化 优 能 力 ， 但 以 上 改进 策略 会 增加 算法 的 计算 量 ， 从 而 增加 算 
算法 。 在 GSA 中 ， 每 个 粒子 的 位 置 表示 问题 的 一 个 可 行 解 ， ”法 的 计算 时 间 ; Joshit1 等 人 在 GSA 中 引入 非 线性 变化 的 引 
解 的 优 劣 与 粒子 的 质量 有 关 ， 质 量 最 大 的 粒子 占据 最 好 的 位 。 力 系 数 ， 提 高 了 GSA 的 寻 优 精度 和 搜索 效率 ， 但 是 没有 对 
。 在 解 空间 内 ， 粒 子 受到 万 有 引力 的 作用 而 相互 吸引 ， 质 ” GSA 的 收 和 伊 性 能 作出 明显 改 ; 
量 小 的 粒子 向 质量 最 大 的 粒子 移动 ， 以 此 完成 算法 的 迭代 寻 综 上 可 知 ， 目 前 已 有 许多 GSA 的 改进 和 应 用 ， 但 GSA 
优 。 基 本 GSA 与 常见 的 生物 启发 式 算法 一 样 , 都 存在 着 易 早 ” 的 全 局 探索 能 力 和 局 部 寻 优 能 力 仍 有 较 大 的 提升 空间 。 为 了 
熟 收 化 、 局 部 搜索 能 力 差 、 缺 少 有 效 加 速 机 制 等 问题 ， 仍 需 。 更 有 效 的 改善 GSA 易 早熟 收敛 的 缺点 ， 增 强 GSA 的 寻 优 性 
对 其 加 以 改进 , 进一步 提高 其 寻 优 性 能 外,。 针对 以 上 问题 , 国 能， 本文 提出 基于 改进 自 适 应 黑洞 机 制 的 引力 搜索 算法 
内 外 学 者 进行 了 大 量 研究 ， 提 出 了 诸多 改进 GSA: 刘 紫 阳 喇 (improved adaptive black hole gravity search algorithm, 
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IABHGSA)。 通 过 改进 Tent 映射 生成 初始 种 群 ， 使 得 种 群 分 
布 随机 、 均 匀 、 遍 历 ， 增 强 算法 的 全 局 勘探 能 力 ， 引入 基于 
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有 二 2E, 0<E<0.5 10 
[20-E) 0.5<E<l (10) 


改进 自 适 应 策略 的 黑洞 机 制 ， 使 得 粒子 可 以 根据 自身 进化 情 
况 更 为 合理 的 选择 位 置 更 新 策略 ， 有 效 降低 算法 陷入 局 部 最 
优 的 可 能 性 ; 引入 基于 学 习 思 想 的 最 优 与 最 差 粒 子 更 新 策略 ， 
增强 算法 的 逃离 局 部 最 优 的 能 力 和 寻 优 速度 ;引入 群体 迁徙 ， 
为 算法 提供 有 效 的 加 速 机 制 ， 提 高 算法 的 收敛 速度 。 最 后 对 
基准 测试 函数 进行 仿真 实验 , 结果 表明 ,， 相 较 于 基本 GSA 和 
其 他 相关 算法 ，IABHGSA 在 寻 优 精度 、 收 敛 速度 及 稳定 性 
上 都 得 到 较 大 的 提升 。 
1 “引力 搜索 算法 

假设 由 VN 个 粒子 X 构成 种 群 , 在 4d 维 引 力 系统 中 ， 定义 
粒子 7 的 位 置 X;=(Xi,X?,.…Xx,.…X”") ， 其 中 i=1,2,.…N 。 根 据 万 
引力 定律 ， 在 a 维 空间 中 ， 第 1 次 迭代 时 粒子 i、 
引力 办 (?) 可 以 表示 为 


FY (1)=G() 


j 之 间 的 


Mi(t)xM,; ; 2 
xy CO) -x¢ (0) 上 


其 中 : G(7) 为 第 1 次 迭代 时 的 引力 系数 ;Mi(n) 与 MD 分 别 为 
粒子 i、j 在 第 1 次 迭代 时 的 惯性 质量 ， 名 (7) 为 第 1 次 迭代 时 


粒子 i、j 之 间 的 欧 氏 距离 ， 为 很 小 的 常量 。 
引力 系数 G(?) 和 欧 氏 距离 乌 () 的 具体 计算 公式 为 
G0)=G,xel 全 (2) 
R, (1)=|X (7), x, (oO) (3) 
其 中 : @ 为 引力 系数 初始 值 ;， 2 为 引力 系数 衰减 因子 ; !、T 
分 别 为 算法 当前 迭代 次 数 和 算法 最 大 办 代 次 数 ， XX(n)、X,(n) 


分 别 为 粒子 i、j 在 第 1 次 迭代 时 的 位 置 。 
粒子 1 的 惯性 质量 M,(n) 根据 适应 度 函 数 定义 为 
Fin, (=min (it, (1) 


mi(?)= max (fis; (1))—min( fiz, (7)) 
mi (1) 
Mi(1)= = 
WF m0 四 


其 中 :应 (0 为 第 1 次 迭代 时 粒子 i 的 适应 度 值 ， max(fis(n))、 
min(fii(n)) 分 别 为 第 1 次 迭代 中 粒子 i 的 最 差 、 最 优 适 应 度 值 。 
为 了 表示 算法 的 随机 性 ， 作 用 在 4 维 上 粒子 i 所 受 合 
(1) 可 以 表示 为 
F(t)= 3 rand; x Fs (1) (6) 


其 中 : randj 为 [0,1] 之 间 的 随机 数 。 在 GSA 的 每 次 迭代 中 ， 粒 
子 的 速度 、 位 置 和 加 速度 都 会 更 新 : 


ve (t+1)=rand, x vy (1)+af (1) (7) 

Xt (t+1)=x¢ (1)+ve (t+1) (8) 
Pr ta) 

af (1)= MD) (9) 


其 中 ， YG) 、 闻 (0 、of 0 分别 为 粒子 i 在 第 ! 代 时 第 4d 维 
间 的 速度 、 位 置 及 加 速度 。 
2 ”改进 引力 搜索 算法 
2.1 改进 Tent 映射 生成 初始 种 群 

基本 GSA 的 初始 种 群 是 随机 生成 的 ， 可 能 导致 个 体 分 
布 不 均 、 种 群 多 样 性 降低 ， 影 响 算法 后 续 的 迭代 寻 优 。 混 沌 
映射 因 其 随机 性 、 人 遍历 性 等 特点 ， 近 年 来 已 被 应 用 于 元 启发 
式 算 法 中 来 提高 种 群 多 样 性 023。 本 文选 用 Tent 映射 来 初始 化 


鼎 


其 中 : 五 为 Tent 映射 。 

分 析 Tent 映射 的 迭代 过 程 能 够 发 现 映 射 中 存在 不 稳定 
周期 点 ， 同 时 由 于 计算 机 的 运行 特性 ， 使 得 Tent 映射 在 经 历 
一 定 次 数 的 迭代 后 总 会 落 入 不 稳定 周期 点 ， 导 致 Tent 映射 的 
遍历 性 降低 。 为 了 避免 Tent 映射 陷入 不 稳定 周期 点 ， 提 高 
Tent 映射 的 均匀 性 和 遍历 性 ， 本 文 提出 基于 logistic 映射 改 
进 的 Tent 映射 ， 利 用 logistic 映射 的 扰动 效果 消除 Tent 映射 


的 不 稳定 周期 点 ， 使 得 Tent 映射 的 分 布 更 人 遍历、 均匀。 本 文 
提出 的 改进 Tent 映射 的 表达 式 为 
Lu=4L(1-L) 0<L<1 (11) 
min DE rand whl 0<E<0.5 
En = ( ) (12) 


min(201—B)+ 南 xrandx "1] QS5<E<l1 


其 中 : 工 为 logistic 映射 ， 5 为 改进 Tent 映射 。 

图 1、2 分别 为 Tent 映射 分 布 直方 图 和 本 文 所 提 改 进 Tent 
央 射 分 布 直方 图 。 从 图 1 中 可 以 看 出 Tent 映射 在 [0,0.1] 范 围 
内 的 取 值 概率 远 高 于 其 他 各 有 段 的 取 值 概率 , 这 是 因为 Tent 映 
射 在 迭代 中 陶 入 不 稳定 周期 点 ， 使 得 后 续 的 映射 值 都 为 0， 
严重 影响 映射 的 均匀 性 和 遍历 性 。 观 察 图 2， 改进 Tent 映射 


在 各 段 的 取 值 概率 相对 均匀 , 表明 本 文 提 出 的 改进 Tent 映射 
成 功 消除 了 使 得 Tent 映射 取 值 为 0 的 点 ， 改 善 了 Tent 映射 
分 布 不 均匀 的 现象 。 
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图 1 Tent 映射 分 布 直 方 图 


Fig.1 Tent map distribution histogram 


number of iimes 


value range 


图 2 改进 Tent 映射 分 布 直 方 图 


Fig.2 Improved Tent mapping distribution histogram 
综 上 ， 本 文 使 用 改进 Tent 映射 生成 初始 种 群 ， 利 用 改进 
Tent 映射 产生 均匀 分 布 的 混沌 序列 ， 减 少 初 始 种 群 分 布 不 均 
匀 对 算法 优化 效果 的 影响 ， 增 强 算法 的 全 局 探索 能 力 。 基 于 
改进 Tent 映射 的 种 群 初始 化 的 表达 式 为 


Xf = Xn 十 (Xu — Xin XE (13) 


种 群 ，Tent 映射 具有 结构 简单 、 遍 历 性 好 、 随 机 性 好 、 混 沌 

均匀 性 和 和 迭代 速度 优 于 logistic 映射 的 特点 03]。 将 Tent 映射 

用 于 GSA 的 种 群 初始 化 , 可 以 增加 初始 种 群 的 多 样 性 ,避免 

算法 陷入 局 部 最 优 ，GSA 可 以 获得 更 好 的 寻 优 效果 。 
Tent 映射 的 表达 式 为 


其 中 :对 为 第 i 个 粒子 在 第 a 维 中 的 位 置 ， a、 xm 分 别 为 

粒子 位 置 的 上 、 下 限 。 

2.2 改进 自 适 应 黑洞 机 制 
zhang 等 人 于 2008 年 首次 提出 随机 黑洞 粒子 群 优化 算法 

(14l(random black hole particle swarm optimization, RBHPSO), 
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RBHPSO 在 PSO 中 引入 黑洞 机 于 
最 优 附 
过 程 。 
局 最 优 gbest 为 球 心 ， 
， 黑 洞 的 吸引 概率 由 [0,1] 之 
如 果 了 <P ， 表 明 粒 子 被 呈 

友之 ， 表 
同 机 制 改进 后 上 
次 ( 人 + 如 (t+1) 了 >P 
gbest +2R(n 
最 优 : R 为 随机 
P 为 黑洞 吸引 概率 。 
必 敛 进程 、 提 
归 洞 机 制 的 GSA(Random 
kt 有 较 强 的 全 


被 吸 纪 
对 于 随机 数 了 ， 
位 置 更 新 公 
用 原 位 置 更 新 公 


优点 ， 所 以 基于 
lgorithm, RBHGSA) 
民 可 能 导致 粒子 过 早 到 达 当 前 


search a 


洞 机 人 


N 式 更 新 粒子 位 置 ; 
。 经 过 黑 让 


wr- 


其 中 : gbest 为 全 局 
之 间 的 随机 数 ; 


许 文 俊 ， 


人 


: 基于 改进 自 适应 黑 
小 使 得 粒子 快速 移动 至 全 后 


近 ， 提 高 了 算法 的 收敛 速度 ， 明 显 加 速 了 算法 的 迭 
黑洞 机 制 假设 在 解 空间 内 存在 一 个 黑洞 ， 该 黑洞 以 
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Xow =g8best + gbest((1—S)xC(0,1)+SxG(0,1)) (19) 


R 为 半径 。 粒 子 进入 黑洞 有 一 定 的 概率 


其 中 : 3 为 自 适 


[gbest f (gbest) < f(x) 
Xbest = | Ka £ (gbest) Sr 大 (Xi ) 


应 权 重 系 数 ; Mew 


车 
党 名 
长 
昌 


间 的 常数 阔 值 P 来 表示 。 
姑 洞 吸引 ， 采 用 新 的 
明 粒 子 逃 离 黑洞 ， 则 使 


的 位 置 更 新 公式 为 


粒子 陷入 


本 文 提出 


其 中 : 


局 部 最 优 的 可 负 
熟 收敛 的 缺点 .为 了 改善 黑洞 机 制 
适应 黑洞 机 制 ， 


ss fit (1)— 应 (1 
| f 


f =max (fis, (1)— fi(t),1) 


-1) Y<P (14) 
黑洞 的 半径 ; Y 、7i 为 [0,1] 


高 算法 收敛 速度 的 


black hole gravity 
局 探索 能 力 。 但 黑 


局 最 优 附 近 ， 增 加 了 
无 法 克服 GSA 易 早 


E 性 ，RBHGSA 
且 易 使 算 


好 (人 + 如 (t+ 


法 早熟 收敛 的 问题 ， 
体 表示 为 


(15) 
(16) 


?<c (17) 


gbest+(2xrand—1)x(gbest—randxx#(t)) ?>c 


名 为 第 ! 代 粒子 的 适应 


1 次 迭代 时 适应 度 值 的 平均 值 ; 


应 度 方差 ; 


万 O] 为 所 有 粒子 在 第 
< 为 适应 度 方差 阔 值 。 
洞 机 制 引 入 适应 度 方 差 吕 取代 黑洞 吸引 


置 ; C(0.1) 为 基 


布 的 变异 算 子 ，f (gbesi)、 


的 最 优 粒 子 的 适应 度 值 ; 


Xbest 为 


(20) 


为 经 过 学 习 后 的 全 局 最 优 位 


于 柯 西 分 布 的 变异 算 子 ，G(0,1) 为 基于 高 斯 分 
了 (wor) 分 别 为 最 优 粒 子 和 经 过 学 习 后 
经 过 选择 后 的 全 局 最 优 位 


基于 自我 学 习 的 最 优 粒子 


更 新 策略 通过 式 (19) 对 全 局 最 


优 进行 自我 学 习 。 在 迭代 初期 ， 


较 大 。 通 过 柯 


尺 后 期 ,1 值 较 大 ， 此 时 高 


小 的 变异 尺度 如 


! 值 较 小 , 此 时 柯 西 算 子 作用 


5 算 子 较 大 的 变异 尺度 对 全 局 最 优 进行 扰动 ， 
避免 算法 落 入 局 部 最 优 ， 增 强 了 算法 的 全 局 探索 能 力 ， 在 先 
斯 算 子 作用 较 大 。 通 过 高 斯 算 子 较 


局 部 范围 内 进行 搜索 ， 算 法 的 收敛 精度 得 到 


提高 ， 在 迭代 中 
好 的 平衡 算法 的 
(20) 选 取 学 习 前 、 


通过 上 述 改进 ， 可 
高 了 算法 的 局 前 


昌 ， 柯 西 算 子 和 高 斯 算 子 出 力 均 多 ， 能 够 更 


后 中 的 较 优 


全 局 探索 能 力 和 局 部 寻 优 能 力 。 最 后 根据 式 
解 作为 新 的 全 局 最 优 , 全 局 最 优 
的 适应 度 值 水 平 得 到 提高 , 以 达到 最 优 粒子 自我 学 习 的 目的 。 
以 有 效 降 低 GSA 陷入 局 
b 寻 优 能 力 ， 算 法 的 收敛 精度 得 到 提高 。 


部 最 优 的 可 能 性 ， 


2.3.2 最 差 粒 子 向 最 优 粒子 的 学 习 


在 TLBO 的 学 习 阶 段 ， 所 有 学 生 


都 会 随机 选择 一 名 同学 


相互 学 习 ， 以 提高 


高 自身 学 习 成 绩 。 


受到 TLBO 中 学 生 相 互 学 


习 的 启发 ， 本 文 提出 了 最 差 粒 子 向 最 优 粒子 的 学 习 策 略 ， 对 


王公 


所 


全 局 最 差 粒子 的 学 习 ， 

Xen = Xuort + (gbest — Xo DLevy(B) CD 
天 op pa de, 
‘Worst We f (Xuom )>f (% 1) 
其 中 :x 为 经 过 


体 表示 为 


(22) 


学 习 后 的 最 差 粒子 的 位 置 ， xwon 为 最 差 粒 


概率 P ,通过 当代 粒子 的 进化 情况 选择 合适 的 位 置 更 新 策略 ， 子 的 位 置 ，@ 表示 为 点 对 点 乘法 ; xx 为 经 过 选择 后 的 全 局 

实现 粒子 位 置 策略 的 自 适 应 选择 。 同 时 改进 原 黑 洞 位 置 最 差 位 置 ， /xs) 、f (Wom) 分 别 为 最 差 粒 子 和 经 过 学 习 后 的 

更 新 公式 ， 以 更 为 随机 的 步 长 向 全 局 最 优 移动 ， 降 最 差 粒 子 的 适应 度 值 ; Levy(Z) 表示 服从 参数 为 8 的 Levy 分 布 

低 粒 子 陷入 局 部 最 优 的 可 能 性 。 当 宁 >c 时 , 表明 粒子 较为 分 散 的 随机 路 径 。 

粒子 陷入 局 部 最 优 的 可 能 性 小 ， 此 时 通过 改进 黑洞 位 置 更 新 策 Levy(O) 表示 如 下 : 

有 信人 局 最 人 也 近 ， 算 流 的 收 信和 到 提 训 ， 当 上 ee 

时 ， 5 集中 ， 粒 子 已 移动 至 全 局 最 优 附近 ， 此 时 通 bp 

过 原 位 置 更 新 策 名 在 局 部 区 域 仔 细 寻 优 ， 寻 优 精 度 得 到 提高 u、v 服 从 正 态 分 布 ， 分 别 表示 为 

综 上 ， 相 较 于 基本 黑洞 机 制 ， 改 进 自 适 应 黑洞 机 制 使 得 粒子 的 u~N(0,02); v~N(0,0?) (24) 

位 置 、 有 效 ， 算 法 可 以 得 到 更 好 的 优化 效果 。 

2.3 关于 学 习 思 息 的 最 优 粒 子 与 最 差 粒 子 更 新 策略 sw 
学 法 [l(teaching-learning-based optimization "(Se Bx27 ee 人 


algorithm， TLBO) 是 Rao 等 人 于 2012 年 提出 的 群 智能 优化 算 


法 ，TLBO 通过 老师 问 学 和 


E 教 学 和 学 生 与 学 生 之 间 相 互 学 习 


两 个 阶段 完成 寻 优 。 本 文 受 到 TLBO 中 教学 与 学 习 的 启发 ， 
融入 到 GSA 的 粒子 更 新 策略 中 ， 


将 TLBO 中 的 教学 与 学 习 

提出 了 基于 学 习 思 想 的 最 

最 优 粒子 的 自我 学 习 和 最 差 粒 子 向 最 优 粒子 的 学 习 ， 
GSA 的 寻 优 精度 和 收敛 速度 。 


2.3.1 基于 自我 学 习 的 最 优 粒子 更 新 

在 TBLO 的 教学 阶段 ， 每 位 学 生根 据 自身 成 绩 向 老师 进 
但 是 随 着 迭代 次 数 的 增 
导致 种 群 多 样 性 迅速 
能 力 。 为 了 提高 算法 的 寻 优 性 
电 的 能 力 ,从 而 得 到 更 好 的 优化 就 果 。 
于 TLBO 中 的 老师 ， 
分 析 ， 本 文 提出 了 基于 自我 学 习 的 最 


其 中 : 工 为 伽 马 函数 。 


最 差 粒 子 向 最 优 粒子 的 学 习 策略 通过 式 (21) 对 全 局 最 差 


粒子 进行 学 习 。 


利用 


Levy 飞行 更 新 全 局 最 差 粒 子 的 位 置 , 加 


和 证 全 


速 全 局 最 差 粒 子 向 全 局 最 优 粒子 移动 ， 减 少 了 算法 的 迭代 时 


行 差异 性 学 习 ， 从 而 
加 ， 老 师 与 学 生 之 间 的 差 明 


提高 全 班 成 绩 。 


丧失 ， 降 低 了 算法 的 局 部 勘 折 


能 ,老师 也 需要 提 


在 GSA 中 ， 最 优 粒 
子 向 其 移动 。 基 于 以 上 


子 相 当 


粒子 更 新 策略 ， 


体 表 示 为 


ssin(£x4) 
T 2 


粒子 与 最 差 粒 子 更 新 策略 ， 通 过 
提 


间 ， 提 高 了 算法 的 收敛 速度 。 同 时 ，Levy 飞行 的 随机 步 长 可 


以 产生 更 随机 的 搜索 过 程 ， 有 利于 保持 种 群 多 样 性 ， 避 免 算 


法 陷入 局 部 最 优 。 


吸引 其 他 粒 


应 度 值 ， 选 择 较 优 
优 粒 子 的 学 习 。 
局 探索 能 力 。 
2.4 群体 迁徙 
在 GSA 寻 优 过 程 中 ， 由 于 惯性 质量 的 累积 效应 Q9， 
由 于 GSA 缺少 有 效 的 加 速 机 4 


强 了 算法 的 全 


通过 上 述 改进 ， 


通过 式 (22) 比 较 学 习 前 、 后 的 最 差 粒 子 的 适 
解 作为 新 的 全 局 最 差 ， 完 成 最 差 粒 子 向 最 


提高 了 算法 的 收敛 速度 ， 增 


前 ， 导 致 粒子 的 速度 随 着 


次 数 的 增加 而 迅速 减 慢 , 削弱 了 GSA 的 全 局 勘探 能 力 , 使 得 


GSA 不 能 有 效 找到 问 


度 ， 提 升 算法 的 划 
群 向 全 局 最 优 收敛 ， 


题 的 最 优 解 。 为 了 增强 GSA 的 收敛 速 
探 性 能 ， 本 文 提 出 群体 迁徙 机 制 ， 加 速 种 
具体 表示 为 
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D2 -x gbesr+ (1+6)x G+D h) 根据 式 (19)(21) 对 最 优 粒 子 和 最 差 粒 子 进行 学 习 ， 并 
2+ rand 按照 式 (20)(22) 更 新 全 局 最 优 和 全 局 最 差 。 

5-rond(-05.03) CI) i) 判断 当前 运算 是 否 满足 停止 条 件 ， 如 果 满 足 ， 输 出 最 

其 中 ,px (+D 为 经 过 群体 迁徙 后 的 粒子 的 位 置 ; 5 为 调控 因子 。 。 优 结果 ， 终 止 运算 ， 若 不 满足 ， 则 跳 转 到 c) 继 续 运算 

群体 迁徙 通过 全 局 最 优 解 和 当前 位 置 引 导 种 群 整体 向 全 ”3.2 时 间 复杂 度 分 析 

最 优 靠近 , 增加 了 算法 的 收敛 速度 , 有 效 改善 GSA 在 途 代 设 算法 种 群 规模 为 V ， 维 度 为 4 ， 最 大 迭代 次 数 为 T 。 
后 期 粒子 移动 速度 慢 的 问题 ， 为 算法 提供 了 一 种 有 效 的 加 速 。 根据 时 间 复 杂 度 符号 o 的 运算 规则 , 对 于 GSA, GSA 初始 化 
收敛 机 制 。 参数 的 时 间 复杂 度 为 00) ， 计 算 GSA 中 各 个 粒子 的 适应 度 


值 的 时 间 复 杂 度 为 O(N)，GSA 完成 了 次 寻 优 的 时 间 复 杂 度 


人 为 O(NxdxT) 。 综 上 ，GSA 总 的 时 间 复杂 度 为 


3.1 算法 流程 O(D+O(N)+O(NxadxT)=O(NxdxT) (28) 
IABHGSA 的 流程 如 图 3 所 示 。 对 于 IABHGSA， 增 加 改进 Tent 映射 生成 初始 种 群 的 时 
间 复 杂 度 为 O(Nxd) ; 增加 改进 自 适 应 黑洞 机 制 的 时 间 复 杂 
ET 度 为 CINWxdxT); 增加 基于 学 习 思 想 的 最 优 与 最 差 粒 子 更 新 
对 种 群 初始 化 策略 的 时 间 复 杂 度 为 0(dxT) ; 增加 群体 迁徙 机 制 的 时 间 复 杂 

根据 适应 度 函 数 计算 所 有 粒子 的 度 为 O(NxdxT)， 则 新 增 改进 策略 的 总 时 间 复 杂 度 为 
站 记 度 信 ” 沦 出 初代 侈 局 昭 估 O(Nxd)+0(NxdxT)+0(dxT)+0(NxdxT) 人 

按照 公式 (3)、(2)、(6)、(9)、(7) 更 新 粒子 的 惯 =O(NxdxT) 

性 质量 、 引 力 系 数 、 所 受 合力 、 加 速度 和 速度 GSA 的 时 间 复 杂 度 为 O(NxdxT) ， 本 文 所 提 改 进 策略 的 
按照 公式 (15) 计 算 当 代 粒 于 适 总 时 间 复 杂 度 为 O(NxdxT) ,所 以 IABHGSA 的 时 间 复 杂 度 为 


应 度 方 差 


O(NxdxT) 。 和 本 文 改进 策略 有 相似 之 处 的 ITC-GSA07 的 时 
间 复 杂 度 为 O(NxdxT) 。 综 上 表明 ， 相 较 于 GSA 和 目前 较 新 
a 0 的 ITC-GSA，IABHGSA 都 没有 增加 计算 负担 。 


粒子 位 置 
4 ”实验 仿真 与 分 析 


为 了 验证 IABHGSA 的 优越 性 能 ， 本 文 设置 五 个 仿真 实 
验 : 第 一 个 实验 将 IABHGSA 与 4 种 近 五 年 内 的 改进 GSA 进 
和 |。 行 比 较 ， 第 二 个 实验 将 IABHGSA 与 种 启发 式 算法 进行 比 

较 ; 第 三 个 实验 将 IABHGSA 与 2 种 改进 RBHGSA 进行 比 
较 ; 第 四 个 实验 将 4 种 结合 单 策 略 的 GSA 与 基本 GSA 进行 
比较 ， 第 五 个 实验 验证 GSA、4 种 改进 GSA、2 种 改进 
RBHGSA 以 及 IABHGSA 在 高 维 测试 函数 上 的 性 能 
4.1 仿真 实验 环境 
图 3 IABHGSA 流程 本 文 在 Intel®@Core™i7-5500U CPU@2.40GHz ， 内 存 8.00GB ， 
Fig.3 Flow chart of IABHGSA Windows10 操作 系统 和 MATLAB 2016a 的 环境 下 对 文中 所 提 所 
= 算法 具体 步骤 如 下 : a) 种 群 初始 化 : 根据 式 (13) 对 种 群 有 算法 进行 仿真 实验 。 
@) 进行 初始 化 : b) 根据 适应 度 函 数 计算 所 有 粒子 的 适应 度 值 ， .2 基准 测试 函数 
选择 出 全 局 最 优 gbest ; c) 按照 式 (5)(2)(6)(9)(7) 更 新 惯性 质量 实验 选取 8 个 不 同类 型 的 基准 测试 函数 ， 包括 4 个 单 峰 
M 、 引 力 系 数 G 、 粒 子 所 受 合力 Ff 、 粒 子 的 加 速度 a 和 速度 "; 基准 测试 函数 、2 个 多 峰 基 准 测 试 函 数 和 2 个 低 维 多 峰 基 准 
d) 按照 式 (15) 计 算 当 代 粒 子 适应 度 方差 避 ; e) 根据 有 F 的 值 选 测试 函数 。 通 过 不 同类 型 的 基准 测试 函数 全 面 验证 
择 式 (8) 或 者 式 (17) 更 新 粒子 位 置 x* ; 了 按照 式 (26) 对 粒子 进行 IABHGSA 的 优越 性 和 本 文 所 提 改 进 策略 的 有 效 性 。 各 个 基 
群体 迁徙 ;，g) 根据 适应 度 函 数 更 新 所 有 粒子 的 适应 度 值 。 准 函 数 的 相关 信息 如 表 1 所 示 。 
表 1 基准 测试 函数 
Tab. 1 Benchmark test functions 


T=T+1 


1 
按照 公式 (26) 进 行 群体 迁徙 


更 新 所 有 粒子 的 适应 度 值 


输出 最 优 结 果 


小 


类 型 基准 测试 函数 范围 维度 ”理论 最 小 值 
帮 =2 好 [-100,100] ”30 0 
,=5" x|+TT x -10,10 30 0 
六 交 f=2 ,5+ I Tx [ ] 
A=2, (2 %) [-100,100] ”30 0 
fi=max {xl,1<i<n} [-100,100] ”30 0 
下 n 六 
-0 池 -Teos( 节 j+1 [-600,600] 30 0 
多 峰 
态 --20o| - E 到 J-en(LE cms) 20+e [-32,32] 30 0 
n 


11 XM(b? +bs) | 
= a A -5,5 4 0.0003 
> b? +bxst+xa | 
维 多 峰 
低 维 多 峰 fi=[1+(m + +1) (19-14% +32 -14x +6%% +3 恕 ) |x 
22.2 2 3 
[30+(25 -3%) x(18-32n +12 好 +48 36 二 27 好 )] 上 


a! 合 { 
hinaz 办 IV 


录用 定稿 许 文 俊 ， 等 : 基于 改进 自 适 应 黑洞 机 制 的 引力 搜索 算法 第 39 卷 第 10 期 


4.3 算法 参数 设置 于 自我 学 习 的 最 优 粒子 更 新 策略 逃离 局 部 最 优 ， 有 效 提高 
本 文选 取 了 4 种 近 五 年 内 的 改进 GSA:CAGSA-PSO081、 算法 的 寻 优 精度 。 


[ou 


SCGSAUH9、IGSAD0、t-IGSAP1H，5 种 启发 式 算 法 :PSOP21、 表 3 算法 寻 优 结果 对 比 
FPA[3] 、BOAPRY、MAPR]、GWORSI 以 及 2 种 改进 Tab.3 Comparison of algorithm optimization results 
RBHGSA:ABHGSAP?"、RBHPSO-GSA[?5] 与 IABHGSA 进行 函数 算法 mean best std 
对 比 实验 。 为 了 使 实验 结果 可 信 、 客 观 ， 设 本 文 所 有 实验 算 GSA 6.82E-08 4.72E-08 1.32E-08 
法 的 种 群 规模 N =30，, 算法 最 大 迭代 次 数 T=1000, 算法 其 他 CAGSA-PSO 2.36E-30 4.21E-31 2.45E-30 
主要 参数 见 表 2。 SCGSA 1.62E-10 2.16E-11 2.84E-10 
表 2 算法 的 主要 参数 . IGSA 2.43E-09 1.07E-09 7.25E-10 
Tab.2 The main parameters of the algorithm t-IGSA 1.09E-18 1.36E-53 5.94E-18 
算法 主要 参数 IABHGSA 2.01E-54 4.76E-71 1.07E-53 
GSA Co =50 ,8g=10 GSA 9.42E-04 7.75E-04 1.23E-04 
CAGSA- y=100,7n=0.1,@ =09 ,wm=04,ci=15,c=0.5 ,其 余 CAGSA-PSO 3.18E-06 1.73E-15 1.32E-05 
PSO 参数 同 GSA SCGSA 3.70E-05 2.31E-05 1.19E-05 
SCGSA Au =2, kwin =0 ,其 余 参 数 同 GSA 1 IGSA 1.63E-04 1.19E-04 2.69E-05 
IGSA x 二 0.9 ,mm 二 0.5 , C1 = 二 0.5 ,Cs =1.5 ,其 余 参 数 同 GSA t-IGSA 2.73E-35 1.78E-45 1.06E-34 
GA | 1/T) ,c=1-1/T ,C=1.8xexp(-1/T) ， 其 余 参 IABHGSA 4.82E-45 5.12E-49 6.19E-45 
数 同 GSA GSA 2.38E+02 1.56E+02 5.90E+01 
PSO Quu =09 , Onin =0.4 ,01 = =2 CAGSA-PSO 3.57E+02 7.90E+01 2.98E+02 
FPA P=rand ,og=0.01,4=1.5 SCGSA 1.38E+02 3.98E+01 4.75E+01 
BOA p=0.6,c=0.01,4a=0.1 1 IGSA 1.12E+00 8.49E-06 2.76E+00 
MA m=1,0=15 ,0 =15, p=2 ,d=5, fi=1 ,LeU(-1l) tIGSA 1.23E+01 8.80E-06 1.93E+01 
GWO Gr =2, dnin =0 IABHGSA 1.21E-19 7.49E-22 1.80E-19 
ABHGSA a=1.5,c,=0.5 ,n=2,p=0.1,R=0.01, 其 余 参 数 同 GSA GSA 1.98E+00 1.05E-04 1.73E+00 
RBHPSO- _] 5 w=05,p=01,R-0.01 ,其 余 参 数 同 GSA CAGSA-PSO 5.51E+00 9.65E-01 2.30E+00 
GSA SCGSA 3.94E-05 2.00E-05 2.41E-05 
IABHOSA c=1x105, B=15 ,其 余 参 数 同 GSA 有 IGSA 3.48E+00 4.62E-01 3.14E+00 
4.4 IABHGSA 与 其 他 改进 GSA 的 比较 t-IGSA 6.16E+00 3.74E-06 3.85E+00 
为 了 验证 IABHGSA 的 优越 性 , 本 节选 取 GSA、CAGSA- IABHGSA 742613 429B-18 119E_12 
PSO、SCGSA、IGSA、t-IGSA 和 IABHGSA 对 8 个 基准 测试 GSA 758E100 346EH00 335EH00 
函数 进行 仿真 实验 。 为 了 使 实验 结果 更 加 客观 ， 每 种 算法 在 CAGSA-PSO 4.76E+00 1.54E400 2 57E+00 
各 个 基准 测试 函数 上 独立 运行 30 次 。 比 较 各 算法 在 独立 的 SCGSA 443E-03 390E13 758E-03 
30 次 运行 中 的 平均 值 (mean)、 最 优 值 (best) 和 标准 差 (std)， 其 广 GX oatod oOo eanpnol 
te 可 以 反映 算法 的 寻 优 精度 ， best 可 以 反映 算法 的 求 tIGSA 1.15E-03 0.00E+00 -3.53E-03 
质量 , std 可 以 反映 算法 的 稳定 性 ,不 同 算法 寻 优 结果 对 比 ， GA GOOBd0. OO DDL 
见 表 3。 GSA 1.79E-04 1.46E-04 2.07E-05 
pa ne 下 到 天 为 单 峰 函数 ， 大 CAGSA-PSO 1.97E-14 7.99E-15 7.05E-15 
至 三 为 多 峰 函 数 。 从 表 2 中 可 以 看 出 ， 无 论 是 求解 单 峰 函数 SCGSA 7.83E-06 3.20E-06 3.16E-06 
种 改进 GSA 和 TABHOSA 党 比 基 本 - h IGSA 3.82E-05 2.41E-05 8.13E-06 
的 寻 人 od 其 中 IABHGSA 的 寻 站 全 全 t-IGSA 4.20E-15 8.88E-16 9.01E-16 
9 种 算法 中 最 好 的 ， 验证 了 IABHGSA 较 强 I IABHGSA 4.09E-15 8.88E-16 1.08E-16 
力 。 在 fi 至 有 上， 部 分 改进 GSA 并 未 起 到 较 好 的 优化 效果 ， Ee Sa 0 
提升 程度 不 大 。 对 于 ，GSA 的 寻 优 结果 与 理论 最 优 相差 较 区 eo dn i 
大 ， 表 明 GSA 在 优化 卢 时 存在 固有 缺陷 :4 种 改进 GSA 对 
卢 的 优化 结果 也 与 理论 最 优 出现 了 较 大 差距 ， 表 明 4 种 改进 卢 es 0 
GSA 并 未 克服 GSA 在 上 上 存在 的 问题 。 而 本 文 所 提 的 IGSA 5.48E-04 3.07E-04 3.56E-04 
IABHGSA 在 六 上 的 寻 优 精度 高 出 GSA 和 4 种 改进 GSA 约 t-IGSA 5.64E-04 3.07E-04 3.36E-04 
20 个 数量 级 ， 极 大 的 提高 了 GSA 在 上 的 寻 优 性 能 ， 克 服 IABHGSA 3.28E-04 3.40E-04 5.05E-05 
了 基本 GSA 在 上 户 上 的 固有 缺陷 。 在 访 至 太 上 ，6 种 算法 的 寻 GSA 3.00E+00 3.00E+00 3.66E-08 
优 结果 的 平均 值 差距 不 大 ,体现 了 GSA 这 类 算法 在 低 维 8 多 CAGSA-PSO 3.00E+00 3.00E+00 1.54E-15 
条 函数 上 具有 较 好 的 寻 优 性 能 。 训 SCGSA 3.00E+00 3.00E+00 5.74E-14 
综 上 ,4 种 改进 GSA 在 一 定 程度 上 提升 了 寻 优 性 能 , 但 IGSA 3.00E+00 3.00E+00 4.40E-11 
是 在 面 对 GSA 的 固有 缺陷 时 无 法 有 效 改 善 。 本 文 所 提 t-IGSA 3.00E+00 3.00E+00 1.28E-14 
IABHGSA 可 以 有 效 提高 算法 的 寻 优 精度 并 且 能 够 有 效 解决 IABHGSA 3.00E+00 。 3.00E+00 2.13E-03 
GSA 存在 的 问题 。 所 以 ， 从 整体 上 看 ，IABHGSA 具有 较 强 为 了 更 直观 的 展示 不 同 算 法 的 寻 优 过 程 ， 本 节 列 出 以 上 
的 竞争 优势 。 这 是 由 于 本 文 所 提 改 进 策略 的 作用 , IABHGSA 种 算法 在 f 至 上 的 迭代 进化 曲线 ， 如 图 4~11 所 示 。 观 察 


通过 改进 自 适应 黑洞 机 制 选择 合适 的 位 置 更 新 策略 ， 使 得 粒 妈 2~7, 可 以 看 出 IABHGSA 具有 不 同 于 其 他 改进 GSA 的 收 
子 位 置 更 新 更 为 合理 、 有 效 。 当 算法 陷入 局 部 最 优 时 ， 通 过 0 在 上 至 大 上 , IABHGSA 在 迭代 前 期 快速 收敛 ,其 收 
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敛 速度 明显 快 于 其 他 5 种 算法 ， 并 且 从 图 中 也 可 以 更 直观 的 1010 

看 出 IABHGSA 的 收敛 精度 得 到 了 较 大 提升 。 对 于 广 和 不 ， 

6 种 算法 都 以 较 快 的 速度 收敛 至 全 局 最 优 附近 ， 收 和 敛 速度 基 ， ET 

本 相同 。 以 上 表明 IABHGSA 提升 了 算法 的 收敛 速度 ， 站 | 一 一 
IABHGSA 具有 和 较 快 的 收敛 速度 主要 是 因为 最 差 粒子 向 最 优 如 cAGSAPSO ~ 

粒子 学 习 策略 和 群体 迁徙 的 作用 。 通 过 最 差 粒 子 向 最 优 粒子 S| de | 
学 习 策略 ， 使 得 最 差 粒 子 快速 移动 至 最 优 粒子 附近 ， 加 速算 算 一 一 tIGSA 

法 的 收敛 。 同 时 由 于 群体 迁移 的 作用 ， 使 种 群 在 位 置 更 新 后 IABHGSA 


整体 向 全 局 最 优 移动 一 定 的 距离 ， 增 强 了 算法 的 全 局 寻 优 能 10% | 
力 。 从 图 4、6、7、9 可 以 看 出 ，IABHGSA 在 迭代 过 程 中 多 


次 跳出 局 部 最 优 ， 在 跳出 局 部 最 优 后 仍 进行 有 效 的 收 化 ， 该 本 | | | | 
过 程 很 好 的 证 明了 本 文 提出 的 改进 策略 可 以 帮助 算法 逃离 局 Yo 200 40 600 800 1000 
部 最 优 , 验证 了 IABHGSA 良好 的 寻 优 性 能 。 综 上 ,IABHGSA 人 迭代 次 数 
提高 了 GSA 的 收敛 速度 ， 加 强 了 GSA 的 全 局 寻 优 能 力 ， 并 图 6 各 算法 对 的 迭代 进化 曲线 
效 克 服 GSA 易 早 熟 收 敛 的 问题 。 Fig.6 Iterative evolution curve of each algorithm to the 六 
标准 差 可 以 检验 算法 的 稳定 性 。 从 表 3 中 可 以 看 出 ， 对 105 [ 


Wy 


于 至 f;，IABHGSA 运行 结果 的 标准 差 是 6 种 算法 中 最 小 
的 , 高 出 其 他 算法 多 个 数量 级 。 对 于 万, IABHGSA 的 稳定 性 


和 其 他 5 种 算法 基本 持平 。 对 于 hh ，IABHGSA 的 稳定 性 较 10" 

差 ， 表明 本 文 提 出 的 改进 策略 并 未 对 GSA 在 大 上 的 运行 产 天 GSA -~ 

生 明显 的 优化 效果 。 根 据 无 免费 午餐 (no-free-lunch, NFL) 定 理 ， 和 和 | 人 ~ - 
没有 任何 一 种 算法 或 者 改进 策略 能 在 所 有 问题 上 达到 最 优 。 一 一 IGSA 

所 以 ， 在 上 IABHGSA 的 稳定 性 较 差 是 可 以 接受 的 。 基 于 二 一 tIGSA 


IABHGSA 


以 上 分 析 ,， 相 较 于 GSA 和 4 种 改进 GSA, IABHGSA 的 稳定 10” 下 
性 得 到 了 提升 ， 算 法 的 鲁 棒 性 增强 。IABHGSA 较 强 的 稳定 


性 得 益 于 最 优 粒 子 的 自我 学 习 , 使 算法 有 能 力 逃 离 局 部 最 优 ， i 所 | | | | | 
在 每 次 寻 优 中 都 可 以 得 到 精度 更 高 的 解 ， 从 而 使 算法 的 稳定 0 200 400 600 800 1000 
性 和 和 鲁 棒 性 得 到 大 幅度 提升 。 适 代数 
图 7 各 算法 对 访 的 迭代 进化 曲线 
Fig.7 Iterative evolution curve of each algorithm to the fs 
100 mg 
过 
Ea 
启 
鉴 地 
I GSA 颖 
和 部 | 一 CAGSA-PSO 攻 i CAGSA-PSO 一 
一 -一 - SCGSA 二 10-10 ea ~ 
一 一 IGSA 一 一 IGSA 了 
= ”mt-IGSA 到 到 tIGSA 
IABHGSA IABHGSA 
0 200 400 600 800 1000 2 
Pe 10- 1 1 1 | 
迭代 次 数 0 200 400 600 800 1000 
图 4 各 算法 对 fi 的 迭代 进化 曲线 迭代 次 数 
Fig.4 Tterative evolution curve of each algorithm to the 大 图 8 各 算法 对 记 的 迭代 进化 曲线 
Fig.8 Iterative evolution curve of each algorithm to the ff; 
105 T T 
100 ED GSA en ,CAGSA- pso 一 "一 " SCGSA 
一 一 JGSA 到 = tIGSA IABHGSA 
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a 
ER 地 
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图 5 各 算法 对 的 迭代 进化 曲线 图 9 各 算法 对 f 的 迭代 进化 曲线 


Fig.5 Iterative evolution curve of each algorithm to the 万 Fig.9 Iterative evolution curve of each algorithm to the 六 
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100 较 高 的 寻 优 精度 。 群 智能 优化 算法 通过 最 优 个 体 的 引导 完成 
GSA 全 局 寻 优 ，IABHGSA 在 此 思想 的 基础 上 引入 了 对 最 优 个 体 

| 的 自我 学 习 ， 使 最 优 个 体 可 以 根据 自身 进化 情况 有 效 更 新 
身 位 置 ， 避 免 由 于 最 优 个 体 的 错误 引导 使 算法 陷入 局 部 最 优 


中 ， 有 效 增强 了 IABHGSA 的 局 部 寻 优 能 力 ， 显 著 提 高 
IABHGSA 的 寻 优 精度 。 
表 4 算法 寻 优 结果 对 比 


Tab.4 Comparison of algorithm optimization results 


示 函 数值 


目 书 


Sr wi i i i ew 


函数 算法 mean best std 

PSO 4.67E-04 3.23E-05 5.49E-04 

ce FPA 1.54E+03 5.23E+01 2.74E+03 

和 迭代 次 数 BOA 1.23E-13 1.14E-13 5.81E-15 

图 10 各 算法 对 六 的 迭代 进化 曲线 MA 4.10E-14 5.58E-16 6.86E-14 

Fig. 10 Iterative evolution curve of each algorithm to the 万 GWO 1.1SE-58 8.76E-61 2.05E-58 

102 ， ， ， ， IABHGSA 2.01E-54 4.76E-71 1.07E-53 

GSA PSO 2.00E+00 2.24E-04 4.22E+00 

人 和 FPA 4.13E+01 2.52E+01 1.16E+01 

tr 

一 | 一 IGSA BOA 6.54E-11 5.88E-11 3.53E-12 

I MA 2.47E-08 3.13E-11 4.23E-08 

各 10 下 ] GWO 1.33E-34 2.39E-35 1.17E-34 

严 上 IABHGSA 4.82E-45 5.12E-49 6.19E-45 

PSO 3.31E+03 3.95E+02 3.49E+03 

FPA 1.51E+04 6.92E+03 7.25E+03 

| BOA 9.84E-14 8.95E-14 9.23E-15 

‘0 a a , MA 1.46B+03 2.36E+02 9.56E+02 

和 迭代 次 数 GWO 8.50E-14 3.28E-20 2.34E-13 

图 11 各 算法 对 的 迭代 进化 曲线 IABHGSA 1.21E-19 7.49E-22 1.80E-19 

Fig. 11 Iterative evolution curve of each algorithm to the _s PSO 5.78E+00 3.99E+00 1.12E+00 

4.5 ”IABHGSA 与 其 他 启发 式 算法 的 比较 FPA 2.43E+00 1.37E+00 1.37E+00 

第 一 个 实验 通过 GSA、4 种 改进 GSA 与 IABHGSA 的 比 BOA 8.11E-11 7.34E-11 4.59E-12 

较 ， 验 证 了 IABHGSA 的 优越 性 能 。 为 了 i 步 验证 村 MA 3.71E+01 2.35E+01 1.03E+01 

IABHGSA 的 优越 性 , 本 节 采 用 5 种 启发 式 算法 与 IABHGSA GWO 1.58E-14 1.23E-15 1.91E-14 

进行 对 比分 析 ， 通 过 不 同类 型 的 算法 比较 进一步 验证 IABHGSA 7.42E-13 4.29E-18 1.19E-12 

IABHGSA 的 优越 性 。 本 节选 用 的 启发 式 算 法 包括 PSO、FPA、 PSO 8.14E-03 3.13E-04 9.61E-03 

CC BOA、MA 以 及 GWO。 选 取 以 上 5 种 算法 对 8 个 基准 函数 FPA 4.25E+01 1.60E+00 7.25E+01 

三 进行 仿真 实验 ， 并 选取 本 文 表 3 中 IABHGSA 的 实验 结果 与 f BOA 2.92E-15 0.00E+00 4.86E-15 

5 种 启发 式 算 法 进行 对 比 。 为 了 使 实验 结果 更 加 客观 ， 每 种 MA 1.08E-02 I.118:15 9.77E-03 

算法 在 各 个 基准 测试 函数 上 独立 运行 30 次 。 比 较 各 算法 在 GWO 1.40E-03 0.00E+00 4.41E-03 

独立 的 30 次 运行 中 的 平均 值 (mean)、 最 优 值 (best) 和 标准 差 IABHGSA 0.00E+00 0.00E+00 0.00E+00 

(std)， 不 同 算法 寻 优 结果 对 比 ， 见 表 4。 PSO 2.57E-02 2.00E-03 3.76E-02 

PSO 是 一 种 经 典 的 启发 式 算法 ， 通 过 群体 中 个 体 自身 的 FPA 1.58E+01 1.01E+01 3.77E+00 

认 知 和 个 体 之 间 的 信息 交流 完成 寻 优 ，PSO 收敛 性 好 ， 结 构 让 BOA 3.86E-11 1.87E-11 1.65E-11 

简单 ， 需 要 设置 的 参数 少 。FPA 模拟 自然 界 中 显 花 植物 花茶 MA 2.32E-01 5.00E-06 4.86E-01 

传粉 过 程 ， 通 过 自 花 授 粉 、 异 花 授 粉 完成 寻 优 ， 具 有 收敛 速 GWO 1.58E-14 1.15E-14 2.80E-15 

度 快 等 优点 。BOA、MA、GWO 是 三 种 较 新 的 基于 种 群 的 优 IABHGSA 4.09E-15 8.88E-16 1.08E-16 

化 算法 , 其 中 BOA 模拟 蝴蝶 的 砚 食 行为 , 通过 香味 的 强度 完 PSO 2.63E-03 3.07E-04 6.25E-03 

成 全 局 搜索 和 局 部 搜索 ，BOA 较 高 的 收敛 精度 ; MA 模 FPA 8.60E-04 6.78E-04 1.75E-04 

拟 雄 、 雌 螺 赃 的 飞行 行为 和 交配 过 程 ， 通 过 交配 产生 下 一 代 f BOA 3.45E-04 3.13E-04 2.22E-05 

完成 寻 优 ，MA 寻 优 能 力 强 、 收 敛 速度 快 ; GWO 模拟 灰 狼 的 MA 6.61E-03 3.07E-04 9.51E-03 

领导 层级 和 狩猎 行为 以 此 寻 优 最 佳 解决 方案 ，GWO 收敛 速 GWO 2.33E-03 3.07E-04 6.34E-03 

度 快 、 寻 优 效 果 好 。IABHGSA 和 上 述 5 种 启发 式 算法 都 属 IABHGSA 3.28E-04 3.40E-04 5.05E-05 

于 群 智 能 优化 算法 ， 其 共同 缺点 为 易 早 熟 收 敛 、 易 陷入 局 部 PSO 3.00E+00 3.00E+00 5.34E-16 

最 优 。 从 表 4 中 可 以 看 出 ，IABHGSA 的 寻 优 效果 是 6 种 算 FPA 3.00E+00 3.00E+00 4.23E-04 

法 中 最 好 的 。 除 IABHGSA 在 A 和 所 上 的 寻 优 结果 ， BOA 3.00E+00 3.00E+00 8.08E-03 

IABHGSA 在 其 他 函数 上 都 得 到 了 6 种 算法 中 的 最 佳 寻 优 精 。 MA 3.00B100 2 3.00E100 2 4.68E-16 

度 。 对 于 及 、 记 、 户 、 有 6，IABHGSA 的 寻 优 精 度 高 出 其 他 篆 GWO 3.00E+00 3.00E+00 1.19E-05 

法 多 个 数量 级 ， 其 寻 优 性 能 的 优势 明显 。 对 于 fh 和 fh， 虽然 IABHGSA 3.00E+00 3.00E+00 2.13E-03 
IABHGSA 的 寻 优 结果 不 是 6 种 算法 中 最 佳 的 ， 但 仍 保持 着 在 稳定 性 方面 , 在 至 f; 上 ,IABHGSA 的 标准 差 是 6 种 
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算法 中 最 低 的 , 稳定 性 最 好 ; 在 户 上 , IABHGSA 的 稳定 性 略 ”4.6 IABHGSA 与 改进 RBHGSA 的 比较 

差 于 BOA， 但 IABHGSA 的 寻 优 精度 优 于 BOA; 在 大 上 ， 为 了 验证 IABHGSA 在 同类 型 算法 中 的 优越 性 ， 本 节选 
IABHGSA 的 稳定 性 不 及 PSO 和 MA， 但 是 在 其 他 函数 上 ” 取 2 种 改进 RBHGSA 以 及 IABHGSA 对 fi 至 进行 仿真 对 
IABHGSA 的 寻 优 精度 和 稳定 性 都 优 于 PSO 和 MA。 以 上 表 比 。 其中, 2 种 改进 RBHGSA 分 别 为 ABHGSA 和 RBHPSO- 
明 ， 虽 然 IABHGSA 在 部 分 函数 上 的 稳定 性 不 及 其 他 算法 ， GSA。ABHGSA 和 RBHPSO-GSA 分 别 在 f 至 上 独行 运行 
但 是 从 总 体 上 看 ，IABHGSA 在 每 次 寻 优 中 都 可 以 得 到 精度 。 30 次 , 得 到 2 种 算法 运行 结果 的 平均 值 (mean)、 最 优 值 (best) 


较 高 的 结果 ， 有 效 提高 了 IABHGSA 优化 成 功 的 概率 ， 使 ” 和 标准 差 (std)。 将 上 述 实 验 结 果 与 本 文 表 3 中 IABHGSA 的 
IABHGSA 的 稳定 性 更 好 、 重 棒 性 更 强 。 实验 结果 进行 对 比 ， 不 同 算法 的 寻 优 结果 对 比 ， 见 表 5。 
图 12~14 为 6 种 算法 在 卢 、 户 、 太 上 的 迭代 进化 曲线 。 在 表 5 算法 寻 优 结果 对 比 
图 12、13 中 可 以 观察 到 PSO、FPA、BOA、MA 因为 收敛 速度 Tab.5 Comparison of algorithm optimization results 
慢 而 得 到 较 差 的 解 。 在 图 14 可 以 看 出 , 虽然 GWO 的 收敛 速度 函数 算法 mean best std 
很 快 ,但 是 其 陷入 了 局 部 最 优 . 从 图 12~14 可 以 看 出 ,IABHGSA ABHGSA 7.73E-09 3.48E-09 ”2.21E-09 
有 较 快 的 收敛 速度 ， 展 示 出 了 良好 的 收敛 特性 ， 并 且 在 迭代 f RBHPSO-GSA 4.12E-09 2.47E-09 1.12E-09 
后 期 具有 较 高 的 寻 优 精度 。 因 此 ， 与 5 种 启发 式 算 法 相 比 ， IABHGSA 2.01E-54 4.76E-71 1.07E-53 
IABHGSA 的 收敛 性 能 较 好 ， 验 证 了 IABHGSA 的 优越 性 。 ABHGSA 4.14E-04 2.49E-04 1.21E-04 
万 RBHPSO-GSA 2.94E-04 1.97E-04 4.82E-05 
IABHGSA 4.82E-45 5.12E-49 6.19E-45 
100 ABHGSA 9.81E-02 2.94E-03 1.31E-01 
fi RBHPSO-GSA 4.03E+00 2.87E-01 5.68E+00 
加 IABHGSA 1.21E-19 7.49E-22 1.80E-19 
a ABHGSA 1.66E-04 1.03E-04 5.09E-05 
严 六 RBHPSO-GSA 3.04E+00 4.58E-05 2.11E+00 
10% | IABHGSA 7.42E-13 4.29E-18 1.19E-12 
ABHGSA 1.71E-02 3.15E-10 2.09E-02 
[ee 天 RBHPSO-GSA 1.56E+00 4.56E-10 1.27E+00 
Od do bd” Bo ~ odd IABHGSA 0.00E+00 0.00E+00 0.00E+00 
和 迭代 次 数 ABHGSA 6.87E-05 5.23E-05 9.98E-06 
图 12 各 算法 对 fi 的 迭代 进化 曲线 fs RBHPSO-GSA 4.75E-05 3.45E-05 7.40E-06 
Fig. 12 Iterative evolution curve of each algorithm to the 六 IABHGSA 4.09E-15 8.88E-16 1.08E-16 
1020 ; . ABHGSA 6.58E-04 3.07E-04 3.08E-04 
万 RBHPSO-GSA 6.03E-04 3.07E-04 2.84E-04 
本 IABHGSA 3.28E-04 3.40E-04 5.05E-05 
ABHGSA 3.00E+00 3.00E+00 7.49E-10 
昌 及 RBHPSO-GSA 3.00E+00 3.00E+00 7.52E-10 
CC 圈 10 | IABHGSA 3.00E+00 3.00E+00 2.13E-03 
一 严 ABHGSA、RBHPSO-GSA 和 IABHGSA 都 是 在 RBHGSA 
1040 | 的 基础 上 作出 改进 ， 基 于 同类 型 算法 的 改进 算法 之 间 的 对 比 
更 能 体现 改进 算法 的 优越 性 和 所 用 改进 策略 的 有 效 性 。 从 表 
JABNGSA 5 中 可 以 看 出 ， 相 较 于 ABHGSA 和 RBHPSO-GSA 的 寻 优 结 
tc 果 ，IABHGSA 在 所 至 上 的 寻 优 结果 都 是 最 好 的 ， 表 明 
迭代 次 数 IABHGSA 在 同类 型 算法 中 竞争 优势 明显 , 验证 了 IABHGSA 
图 13 各 算法 对 户 的 迭代 进化 曲线 有 较 强 的 寻 优 性 能 。 同 时 从 表 5 中 可 以 看 出 ， 对 于 使 用 基 
Fig. 13 Iterative evolution curve of each algorithm to the 万 本 随机 黑洞 的 ABHGSA 和 RBHPSO-GSA， 其 寻 优 结果 相 较 
于 表 3 中 GSA 的 实验 结果 , 提升 程度 有 限 ,， 大 概 提高 了 1 个 
数量 级 左右 。 这 是 因为 基本 随机 黑洞 策略 只 能 随机 地 在 两 种 
10° 立 置 更 新 策略 中 进行 选择 ， 这 样 会 导致 算法 过 早 抵达 当前 全 
局 最 优 附 近 ， 从 而 早熟 收敛 ， 削 弱 了 算法 的 局 部 勘探 能 力 。 
通 IABHGSA 不 同 于 以 上 两 种 改进 RBHGSA，IABHGSA 对 基 
本 随机 黑洞 策略 进行 改进 ， 根 据 当 代 粒 子 适应 度 方差 选择 合 
全 nd | 适 的 位 置 更 新 策略 ， 使 粒子 位 置 的 更 新 更 符合 粒子 的 迭代 情 
况 ， 避 人 免 因 随机 选择 位 置 更 新 策略 而 导致 的 算法 陷入 局 部 最 
优 ， 从 而 有 效 提 高 了 算法 的 局 部 寻 优 能 力 。 同 时 ， 基 于 自我 
学 习 的 最 优 粒子 更 新 策略 可 以 避免 算法 陷入 局 部 最 优 ， 也 降 
辣 PA | | 低 了 算法 陷入 局 部 最 优 的 可 能 性 。 两 种 改进 策略 使 
0 200 400 600 800 1000 IABHGSA 的 寻 优 精度 得 到 较 大 的 提升 ,验证 了 IABHGSA 在 
和 迭代 次 数 同类 型 算法 中 优异 的 寻 优 性 能 。 在 算法 的 稳定 性 方面 ， 
图 14 各 算法 对 户 的 迭代 进化 曲线 IABHGSA 的 综合 表现 较 好 ， 仅 在 大 上 的 稳定 性 不 如 2 种 改 


Fig. 14 Iterative evolution curve of each algorithm to the 广 进 RBHGSA。 综 上 表明 ，IABHGSA 的 稳定 性 在 超过 85% 的 
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测试 函数 上 均 有 明显 的 优势 ， 算 法 的 寻 优 效果 更 为 稳定 ， 在 从 表 6 中 可 以 看 出 ， 在 至 上，4 种 结合 单 策略 改进 
面 对 各 种 优化 问题 时 的 稳定 性 和 和 鲁 棒 性 更 好 。 的 GSA 的 寻 优 精度 和 稳定 性 都 得 到 了 提升 ， 表 明 结 合 单 策 
4.7 不 同 改进 策略 的 比较 各 改进 的 GSA 的 寻 优 性 能 和 重 棒 性 更 好 。GSA2 和 GSA3 在 
本 节 通 过 GSA 与 4 种 结合 单 策 略 改进 的 GSA 进行 比较 ， fi 至 上 的 寻 优 效 果 较 好 ， 这 验证 了 改进 自 适应 黑洞 机 制 和 
从 而 验证 本 文 所 提 改 进 策 略 的 有 效 性 。 将 结合 改进 Tent 映射 ” 基于 自我 学 习 的 最 优 粒子 更 新 策略 可 以 有 效 改 善 算 法 易 早熟 
生成 初始 种 群 的 GSA 命名 为 GSA1; 结合 改进 自 适 应 黑洞 机 ”收敛 的 缺点 ， 提 升 算法 的 寻 优 精度 。 对 于 方 和 有 ，GSA 和 4 


制 的 GSA 命名 为 GSA2; 结合 基于 学 习 思 想 的 最 优 与 最 差 粒 。 种 结合 单 策略 改进 的 GSA 的 寻 优 精度 基本 相同 。 在 稳定 性 


子 更 新 策略 的 GSA 命名 为 GSA3; 结合 群体 迁徙 的 GSA 命 上，GSA2 在 和 上 表现 不 佳 ， 但 GSA2 在 其 他 函数 上 的 
名 为 GSA4。 将 GSA 与 4 种 结合 单 策略 的 GSA 进行 对 比 实 收敛 精度 和 稳定 性 都 优 于 GSA， 根 据 NFL 定理 ， 由 此 说 明 
验 ， 并 给 出 GSA1、GSA2、GSA3、GSA4 在 fi 至 上 独立 运 GSA2 中 的 改进 策略 仍 是 有 效 的 。 
行 30 次 后 得 到 的 寻 优 结果 的 平均 值 (mean)、 最 优 值 (best) 和 从 图 15~18 中 可 以 看 出 , 4 种 单 策略 改进 的 GSA 的 收敛 
标准 差 (std), 与 本 文 表 3 中 GSA 的 实验 结果 进行 对 比 ， 相 关 ” 速度 优 于 GSA, 表明 本 文 所 提 的 改进 策略 可 以 有 效 提高 算法 
数据 见 表 6。 由 于 篇 幅 原因 ， 但 为 了 更 直观 的 比较 不 同 改进 的 收敛 速度 , 增强 了 GSA 的 全 局 寻 优 能 力 , 对 于 图 17, GSAI1 
策略 对 GSA 收敛 速度 的 影响 ， 本 节 给 出 各 算法 对 hh、 所 、 的 收敛 速度 得 到 提升 , 验证 了 改进 Tent 映射 生成 初始 种 群 可 
廊 、 记 的 迭代 进化 曲线 ， 如 图 15~18 所 示 。 以 增强 算法 的 全 局 探索 能 力 。 观 察 图 17 和 18,GSA2 和 GSA3 
表 6 算法 寻 优 结果 对 比 的 曲线 上 存在 多 个 转折 点 ， 表 明 改 进 自 适 应 黑洞 机 制 和 基于 
Tab.6 Comparison of algorithm optimization results 我 学 习 的 最 优 粒 子 更 新 策略 可 以 有 效 的 帮助 算法 逃离 局 部 
函数 算法 Imean best std 最 优 ， 提 升 算 法 的 寻 优 精度 。 同 时 ，GSA4 在 i、 方 和 万 上 
GSA 6.82E-08 4.72E-08 1.32E-08 收敛 速度 提升 明显 ， 表 明 群 体 迁 徙 策略 可 以 加 速算 法 收敛 ， 
GSA1 7.65E-09 5.19E-09 2.54E-09 为 算法 提供 了 有 效 的 加 速 收敛 机 制 ， 克 服 了 GSA 缺少 加 速 
GSA2 1.10E-17 2.15E-19 2.21E-17 机 制 的 缺点 。 
GSA3 1.10E-41 6.94E-51 3.23E-41 总 结 而 言 , GSA3 的 寻 优 结果 的 有 效 性 显著 , 在 寻 优 精度 、 
GSA4 1.89E-12 1.10E-15 3.71E-12 收敛 速度 、 稳 定性 上 都 有 较 大 提升 。GSA1 的 寻 优 结果 是 4 种 
GSA 9.42E-04 7.75E-04 1.23E-04 单 策 略 改进 的 GSA 中 效果 最 差 的 , 但 其 三 个 评价 指标 都 略 优 于 
GSAI1 4.31E-04 3.10E-04 5.70E-05 GSA。GSA2 在 方 和 大 上 表现 较 差 ， 但 在 其 他 函数 上 的 寻 优 精 
fh GSA2 6.82E-13 1.31E-13 5.98E-13 度 、 收 敛 速度 、 稳 定性 相 较 于 GSA 得 到 了 较 大 提升 。GSA4 在 
GSA3 3.92E-39 5.44E-47 5.92E-39 户 上 收敛 速度 与 GSA 持平 ,但 在 AA、 有 、 万 上 收敛 速度 高 于 
GSA4 8.84E-07 2.05E-07 1.20E-06 GSA， 并 且 GSA4 在 fi 至 上 有 较 高 的 寻 优 精 度 和 稳定 性 。 综 
GSA 2.38E+02 1.56E+02 5.90E+01 上 ， 虽 然 部 分 改进 策略 在 某 些 函数 上 表现 不 佳 ， 但 是 从 总 体 上 
GSAI1 1.95E+01 1.74E+00 1.77E+01 看 是 有 效 的 ， 验 证 了 本 文 所 提 改 进 策略 的 有 效 性 。 
基 GSA2 3.13E-19 8.18E-21 4.34E-19 
GSA3 6.87E+01 3.26E+01 3.38E+01 
GSA4 1.11E-12 3.20E-16 3.08E-12 
GSA 1.98E+00 1.05E-04 1.73E+00 
GSA1 8.30E-01 4.01E-05 2.13E+00 
记 GSA2 2.02E-12 5.40E-14 3.51E-12 全 
GSA3 5.77E-01 3.69E-05 9.51E-01 严 
GSA4 7.52E-07 5.61E-09 8.98E-07 
GSA 7.58E+00 3.46E+00 3.35E+00 
GSA1 2.89E+00 2.27E-10 5.50E+00 
六 GSA2 0.00E+00 0.00E+00 0.00E+00 5 06 Ee Bo — oDd 
GSA3 0.00E+00 0.00E+00 0.00E+00 和 迭代 次 数 
GSA4 3.46E-04 0.00E+00 1.89E-03 图 15 各 算法 对 f 的 迭代 进化 曲线 
GSA 1.79E-04 1.46E-04 2.07E-05 Fig. 15 Iterative evolution curve of each algorithm to the fi 
GSA1 6.74E-05 4.15E-05 1.09E-05 1020 
区 GSA2 1.71E-12 5.06E-14 1.86E-12 
GSA3 4.09E-15 8.88E-16 1.43E-15 
GSA4 2.48E-07 9.81E-09 4.81E-07 
GSA 2.45E-03 8.22E-04 1.27E-03 3 
GSA1 1.88E-03 6.25E-04 2.24E-03 乓 
fh GSA2 1.75B-02 7.41E-04 2.99E-02 严 | a 、、 
GSA3 1.05E-03 5.11E-04 842E-04 | GSA1 \ 
ma en 让 二 
GSA4 3.63E-04 3.08E-04 2.42E-04 GSA 六 
GSA 3.00E+00 3.00E+00 3.66E-08 GSA4 
GSAl 3.00E+00 3.00E+00 3.19E-10 cc 
六 GSA2 2.06E+01 3.01E+00 2.78E+01 选 代 次 数 
GSA3 3.00E+00 3.00E+00 2.70E-10 图 16 各 算法 对 万 的 迭代 进化 曲线 
GSA4 3.00E+00 3.00E+00 3.87E-04 Fig. 16 Iterative evolution curve of each algorithm to the 户 
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于 基本 GSA 和 6 种 改进 GSA。 以 上 表明 IABHGSA 可 以 有 
效 的 处 理 复杂 的 高 维 问题 ， 验 证 了 IABHGSA 的 优越 性 和 本 
文 所 提 改 进 策略 的 有 效 性 。 
表 7 算法 对 高 维 测试 函数 寻 优 结果 对 比 


3 Tab.7 Comparison of optimization results of algorithms on 
要 high-dimensional test functions 
本 函数 算法 mean best std 
GSA 4.06E+04 3.66E+04 3.31E+03 
CAGSA-PSO 1.30E+05 9.71E+04 1.98E+04 
SCGSA 1.85E+04 1.56E+04 2.29E+03 
cc IGSA 1.61E+05 1.13E+05 2.18E+04 
和 迭代 次 数 4 t-IGSA 2.02E+05 1.67E+05 4.05E+04 
图 17 各 算法 对 f; 的 迭代 进化 曲线 ABHGSA 5.72E+04 4.65E+04 8.99E+03 
Fig. 17 Iterative evolution curve of each algorithm to the 态 RBHPSO-GSA 1.40E+05 1.19E+05 2.28E+04 
IABHGSA 1.07E-17 4.21E-20 2.70E-17 
GSA 5.23E+266 1.59E+02 Inf 
CAGSA-PSO 5.23E+266 2.87E+02 Inf 
SCGSA 5.23E+266 2.57E+02 Inf 
IGSA 5.23E+266 3.78E+02 Inf 
ee 问 2” t-IGSA 2.57E+01 4.60E-11 4.25E+01 
严 ABHGSA 1.71E+02 1.53E+02 1.49E+01 
RBHPSO-GSA 2.48E+02 1.95E+02 3.22E+01 
IABHGSA 2.56E-04 9.87E-08 4.99E-04 
GSA 2.84E+05 1.74E+05 1.36E+05 
ce CAGSA-PSO 1.40E+06 7.76E+05 6.78E+05 
迭代 次 数 SCGSA 1.50E+05 1.07E+05 3.11E+04 
图 18 各 算法 对 记 的 迭代 进化 曲线 IGSA 3.59E+06 2.05E+06 1.26E+06 
Fig. 18 Iterative evolution curve of each algorithm to the 万 t-IGSA 3.13E+06 5.46E+05 2.00E+06 
4.8 ”高 维 测试 函数 实验 结果 及 分 析 ABHGSA 1.04E+06 4.61lE+05 7.47E+05 
本 文 所 提 IABHGSA 在 低 维 函数 上 有 较 好 的 寻 优 精度 、 RBHPSO-GSA = 2.18E+06 7.75E+05 1.41E+06 
收敛 精度 和 稳定 性 。 但 是 一 般 的 改进 策略 在 面 对 实际 工程 问 IABHGSA 7.53E-19 4.51E-20 1.00E-18 
题 ， 特 别 是 高 维 测试 函数 时 极 易 失效 。 为 了 进一步 验证 GSA 2.58E+01 2.31E+01 2.03E+00 
IABHGSA 的 优越 性 和 本 文 所 提 改 进 策略 的 有 效 性 ， 本 节选 CAGSA-PSO 4.74E+01 4.46E+01 3.11E+00 
取 GSA、CAGSA-PSO、SCGSA、IGSA 、t-IGSA、ABHGSA、 SCGSA 2.22E+01 2.03E+01 1.71E+00 
RBHPSO-GSA 以 及 IABHGSA 对 fi 至 所 进行 高 维 测试 函数 IGSA 5.99E+01 5.11E+01 4.00E+00 
的 仿真 实验 .高 维 维度 设置 为 500 维 , 关于 其 他 参数 的 设置 ， t-IGSA 5.69E+01 4.83E+01 5.43E+00 
己 在 本 文 4.3 节 中 详细 阐述 。 为 了 使 实验 结果 更 加 客观 ， 每 ABHGSA 4.57E+01 4.20E+01 3.26E+00 
种 算法 在 各 个 高 维 测试 函数 上 独立 运行 30 次 。 比 较 各 算法 RBHPSO-GSA 4.92E+01 4.43E+01 3.93E+00 
在 独立 的 30 次 运行 中 的 平均 值 (mean)、 最 优 值 (best) 和 标准 IABHGSA 7.78E-12 1.73E-12 4.26E-12 
差 (std)， 不 同 算法 寻 优 结果 对 比 ， 见 表 7。 GSA 4.57E+03 4.40E+03 1.15E+02 
从 表 7 中 可 以 看 出 ， 在 高 维 函 数 上 ，IABHGSA 的 收敛 CAGSA-PSO 1.65E+03 1.42E+03 1.66E+02 
精度 依旧 较 高 ， 在 至 上 都 收敛 至 理论 最 优 附近 ， 表 明 SCGSA 4.03E+02 3.53E+02 3.41E+01 
IABHGSA 在 高 维 函 数 上 依然 有 较 好 的 寻 优 性 能 。 对 于 GSA 5 IGSA 2.01E+03 1.31E+03 4.64E+02 
和 6 种 改进 GSA， 其 失效 现象 较为 明显 ,尤其 在 有 上。 是 | t-IGSA 2.14E+03 1.92E+03 1.20E+02 
连续 的 、 平 滑 的 单 模 态 函数 ， 当 自 变量 趋 近 于 无 穷 大 时 ， 函 ABHGSA 1.49E+03 1.14E+03 1.95E+02 
数 会 形成 大 量 局 部 极 值 区 域 ， 易 陷入 局 部 最 优 且 不 易 跳出 ， RBHPSO-GSA 1.81E+03 1.43E+03 2.50E+02 
在 高 维 测试 函数 下 求解 难度 较 大 。GSA、CAGSA-PSO、 IABHGSA 3.33E-16 0.00E+00 6.62E-16 
SCGSA、IGSA 在 户 上 严重 失效 ， 几 乎 未 起 到 优化 效果 ; 虽 GSA 9.16E+00 ”8.71E+00 2.91E-01 
然 tIGSA、ABHGSA 和 RBHPSO-GSA 的 寻 优 结果 优 于 GSA、 CAGSA-PSO 1.35E+01 1.25E+01 7.71E-01 
CAGSA-PSO、SCGSA、IGSA, 但 是 失效 现象 依旧 较为 明显 ， SCGSA 8.51E+00 7.75E+00 5.11E-01 
这 也 证 明了 一 般 的 改进 策略 在 高 维 函 数 上 极 易 失效 。 IGSA 1.54E+01 1.46E+01 4.56E-01 
IABHGSA 在 上 的 收敛 精度 高 出 GSA 和 6 种 改进 GSA 多 6 tIGSA 1.55E+01 1.46E+01 5.47E-01 
个 数量 级 ， 收 敛 精度 远 高 于 GSA 和 6 种 改进 GSA， 这 是 因 ABHGSA 1.29E+01 1.12E+01 7.65E-01 
为 基于 自我 学 习 的 最 优 粒子 更 新 策略 可 以 帮助 算法 逃离 局 部 RBHPSO-GSA 1.32E+01 1.23E+01 5.90E-01 
最 优 ， 增 强 算法 的 局 部 寻 优 能 力 ， 提 高 算法 的 收敛 精度 。 同 IABHGSA 1.12E-11 1.57E-12 6.38E-12 
时 ，IABHGSA 在 fi 至 6 寻 优 结果 的 标准 差 都 低 于 其 他 7 种 5 ”结束 语 
算法 ， 表 明 在 高 维 函数 上 IABHGSA 的 稳定 性 依旧 较 强 、 重 SH 
棒 性 好 。 综 上 , 对 于 高 维 函数 , IABHGSA 的 寻 优 性 能 依旧 强 本 文 提出 了 一 种 基于 改进 自 适 应 黑洞 机 制 的 引力 搜索 算 


许 文 俊 ， 


法 ,使 用 基于 logistic 映射 改进 的 Tent 映射 进行 种 群 初始 化 ， 

改善 算法 的 全 局 勘探 能 力 。 引 入 改进 自 适应 黑洞 机 制 ， 使 粒 
子 按照 和 欠 代 进化 情况 选择 位 置 更 新 策略 ， 使 得 位 置 更 新 更 为 
合理 ， 同 时 提升 了 算法 的 寻 优 精度 。 引 入 基于 学 习 思 想 的 最 
CE 与 最 差 粒 子 更 新 策略 ， 有 效 避 人 免 算法 陷入 局 部 最 优 和 提升 
算法 的 收敛 速度 。 采 用 群体 迁徙 机 制 ， 有 效 增强 了 算法 的 全 
局 探索 能 力 ， 为 算法 提供 了 有 效 的 加 速 机 制 。 最 后 ， 通 过 仿 
验证 了 IABHGSA 的 优越 性 和 本 文 所 提 改 进 策 
由 于 IABHGSA 的 研究 处 于 初始 阶段 ， 接 下 来 的 研究 方向 主 
要 为 拓展 IABHGSA 至 实际 问题 中 的 应 用 。 
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